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We present a series expansion study of spin-S square-lattice Heisenberg antiferromagnets. The nu- 
merical data are in excellent agreement with recent neutron scattering measurements. Our key result 
is that the correlation length £ for S > 1/2 strongly deviates from the exact T — > (renormalizcd 
classical, or RC) scaling prediction for all experimentally and numerically accessible temperatures. 
We note basic trends with S of the experimental and series expansion correlation length data and 
propose a scaling crossover scenario to explain them. 



In recent years much attention has focused on two- 
dimensional (2D) square-lattice quantum Heisenberg an- 
tiferromagnets, described by the Hamiltonian 
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where (ij) denotes summation over all pairs of near- 
est neighbors. In their seminal work [jjj, Chakravarty, 
Halperin, and Nelson (CHN), utilizing a mapping of the 
low-energy spectrum of Eq. (jl|) onto the quantum non- 
linear sigma model (QNLerM), have shown that the low- 
temperature properties of these systems obey renormal- 
ized classical (RC) scaling, where the correlation length 
£ w 0.37 x (c/p s )exp(2irp s /T) for T < p s (here p s is 
the T — > spin stiffness, and c is the spin- wave velocity). 
Subsequently, Hasenfratz and Niedermayer calculated for 
the 2D QNLerM the exact value of the prefactor and the 
leading 0(T/2np s ) correction H: 
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The neutron scattering measurements of £(T) in 
5 = 1/2 layered Heisenberg antiferromagnets such as 
La2Cu04 M and Sr2Cu02Cl2 M reveal a remarkable 
agreement with Eq.(||). At first sight, one would expect 
the RC description to improve as the value of the spin 
increases. If 5 is formally regarded as a continuous vari- 
able, then the Neel order is expected to vanish for some 
5 < 1/2. At the critical point, where p s vanishes, Eq.(||) 
fails and the correlation length is, in fact, inversely pro- 
portional to the temperature; this is the quantum critical 
regime (QC), discussed in Ref. @,| (it was argued ear- 
lier |||| that the 5=1/2 Heisenberg model ([!]) exhibits 
certain signatures of QC behavior for T > 0.5J; we do 
not discuss RC to QC crossover effects in this letter). 
Naively, increasing the value of spin moves the system 



away from this limit so that the RC behavior would be 
more pronounced. 



FIG. 1. Semi-log plot of the series expansion result for the 
correlation length in the 5=1 Heisenberg antiferromagnet 
versus T/J (solid line) together with the experimental data 
for K2NiF4 (solid circles, ||) and La2Ni04 (open circles, j?|]). 
Numerical and experimental results agree with each other, 
and all deviate from the exact RC prediction |B| evaluated 
using known values for c and p s (dashed line). 

However, such an expectation does not hold. Re- 
cently, Greven and co-workers QJ^] have reported a sig- 
nificant discrepancy between the neutron scattering mea- 
surements of the correlation length in the 5=1 systems 
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K2N1F4 and La2Ni04 and the RC prediction. Prelimi- 
nary experiments on the S — 5/2 system Rb2MnF4 reveal 
an even larger discrepancy . As is evident from Fig j|, 
our series expansion result for S = 1 |ll]] is in excel- 
lent agreement with the experimental data in the region 
of overlap, and also deviates substantially from the RC 
prediction (note, that £ is expressed in units of the lattice 
constant a). 

In order to investigate the origin of this deviation, 
we have calculated high-temperature expansions for the 
Fourier transform of the spin-spin correlation func- 
tion (StqSq) for all spins in the range S — 1/2 to 
S = 5/2. The generation of the series expansions 
will be discussed elsewhere. We present the data for 
the T = structure factor Sq — S(Q), where Q = 
(■7r/a,7r/a) is the staggered ordering wavevector, and 
the second moment correlation length £ is defined by 
£ 2 = - (<9 2 S(Q + q)/dq 2 ) I /2Sq. The series are ana- 
lyzed by either performing direct Pade approximation or 
by taking the logarithm of the series first and then calcu- 
lating the Pade approximants. The latter is likely to show 
better convergence if the correlation length increases ex- 
ponentially fast at low temperatures. We restrict our- 
selves to temperature ranges where different methods of 
extrapolation agree within a few per cent. 



The ratio between the calculated £ and the Hasenfratz 
and Niedermayer formula £hn given by Eq. (^) , is plotted 
in Fig.|| for different spins as a function of T/p s . Con- 
trary to the naive expectation that the RC behavior be- 
comes more pronounced as S increases, we find that for 
the range of temperatures probed here £ monotonically 
deviates from £hn as S increases. As noted above, this 
same systematic trend is seen experimentally [pPj7|,^0| . 



FIG. 2. Series expansion results for the correlation length 
plotted as £/£hn versus T/p a . Here £hn is the exact RC 
prediction g], Eq.(^), with c and p a from Ref. |^|]]. The ratio 
£/£hn, which is fairly close to unity for S = 1/2, gradually 
decreases as S increases, suggesting increasing disagreement 
with RC theory for larger S. 



FIG. 3. Semi-log plot of S£/a versus T/JS 2 of 
Sr 2 Cu0 2 Cl 2 (S = 1/2, open circles, [§) and K 2 NiF 4 (S = 1, 
solid squares, [Q) together with the series expansion results 
(1/2 < S < 5/2, represented by the various lines) as well as 
Monte Carlo data (S = 1/2, solid squares, JT^). However 
impressive, this data collapse is inconsistent with Eq.(Q) with 
c and p s from Ref. 

In an attempt to understand these data, we have 
considered various possible scaling scenarios. To high- 
light the expected spin dependence, we note that for a 
square-lattice nearest-neighbor Heisenberg antiferromag- 
net, p s = Z p {S)JS 2 and c = 2 3 / 2 aZ c (S) JS. The quan- 
tum renormalization factors Z p (S) and Z C (S) are known 
from spin- wave theory |§ (for all values of S), T = 
series expansion II (for S = 1/2 and S — 1), and Monte 
Carlo studies @ (for S = 1/2). For S = 1/2 and S = 1 
these different methods yield good agreement. 

Eq.(||) may then be written as 

S£ eZ c ( 2nZ p \ \ T/JS 2 1 

which suggests plotting S£/a versus T/JS 2 to elucidate 
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the dependence on 5. The results so-obtained are shown 
m Fig]| For S = 1/2 we also show Monte Carlo data 
jl2| , which extend to somewhat lower temperatures than 
the series expansion result. Surprisingly, over the range 
of 5£ I a from 1 to at least 50, the data to a good ap- 
proximation fall on the same curve. For the sake of 
clarity, we have omitted experimental data for La2Cu04 
(5 = 1/2), La 2 Ni0 4 (5 = 1), and Rb 2 MnF 4 (5 = 5/2), 
all of which fall on the approximate "scaling" curve of 
Fig. m to within experimental error. 

It was demonstrated in Ref. [Q that Eq.(||) describes 
the correlation length data in absolute units in the 
5 = 1/2 system S^CuC^C^ extremely well. Inter- 
preted naively, Fig.[| would then suggest that RC be- 
havior holds for all S, but with quantum renormalization 
factors Z p (S) and Z C (S) that are nearly S-independent 
for 1/2 < 5 < 5/2 and close to their values at S = 1/2, 
Z p (l/2) w 0.72 and Z c (l/2) « 1.2. However, this is very 
unlikely. Spin-wave theory @] predicts a substantial S- 
dependence: Z p (5/2) « 0.95 and Z c (5/2) « 1.03, which 
is already close to the classical limit (5 — > oo) in which 
Z p , Z c — » 1. Serious errors in these values seem very un- 
likely, given their good agreement for S = 1/2 and 5 = 1 
with those obtained by other methods [p|Jl^] . 

In an attempt to understand these data, it is help- 
ful to recall that at any fixed T/p s RC theory will in- 
evitably fail for sufficiently large values of spin. Indeed, 
a straightforward application of Eq.(|^) to the 5 — > oo 
limit taken at T/JS 2 ~ 1 would predict £ — 0. However, 
the above limit corresponds to the classical Heisenberg 
magnet, where £/a is known to be finite and of order 
unity for T ~ JS 2 . 

One can understand where Eq. (^J) may fail by following 
CHN in their derivation of the leading asymptotic behav- 
ior in the RC regime, but taking into account that 5 may 
be large. CHN have shown that for T <C p s the magnetic 
correlations can be calculated using classical dynamics, 
except that all wavevector integrations should be limited 



to |q 



< „. = 



T/c rather than taken over the whole Bril- 



louin zone. Here the words "classical dynamics" simply 
mean that for |q| < q c , all Bose factors for spin-waves 
can be approximated assuming cq <C T. The key result 
of CHN is to show that such a calculation yields the cor- 
rect spin correlations for the quantum Heisenberg model. 
We now evaluate the cutoff wavevector as 
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This dependence of q c on 5 arises because the spin stiff- 
ness is proportional to the square of 5, but the spin- wave 
velocity is only linear in 5. 

For S > 1 and T ~ p s ~ JS 2 , the cutoff wavevector 
q c ~ S/a it I a is outside the Brillouin zone. Hence, 
the requirement that cq < T, or equivalently q < q c , 
places no further restrictions on the q-integrations which 
are already limited by the Brillouin zone. In this case 



all of the integrals are the same as those of the classical 
Heisenberg magnet, and the classical 5 — ► oo limit is 
recovered. 

The crossover temperature T cr between the RC and 
classical regimes depends on 5, and its order of magni- 
tude can be estimated as the temperature where q c ~ 
a -1 . This yields 



JS, while p s ~ JS 2 
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By substituting T cr into Eq.(g), one concludes that the 
crossover from RC behavior at low temperatures to clas- 
sical behavior at higher temperatures should occur for a 
£cr = £(T cr ) that is larger for larger 5. 



FIG. 4. Semi-log plot of the correlation length ob- 
tained from series expansion plotted versus T/(JS(S + 1)) 
for S = 1/2, 1, and 5/2. For larger spins £/a is close to the 
classical (S — » oo) limit, which provides evidence that classi- 
cal (S — > oo) magnetic behavior holds for JS < T < JS 2 , 
in agreement with our proposed scaling crossover scenario. 
Note, that in most of the temperature range shown, the 
S — > oo model is not in the scaling limit and its corre- 
lation length deviates from the expected T — > behavior 
£/a ~ const x (T/JS 2 ) exp(27r JS 2 /T). 

In order to test this scenario, we plot the correlation 
length as a function of T/{JS{S + 1)) in Fig.|, where 
J5(5 + 1) is the classical (not T = 0) spin stiffness. In 
replacing 5 2 by 5(5 + 1), we follow a purely empirical 
observation that the correlation length at T > JS for 
5 > 3/2 seems to depend on 5 primarily through the 
combination 5(5+ 1) (for T <C JS, £ depends on 5 only 



3 



through p s and c). We find that £ gradually approaches 
the classical result as 5* increases, and the difference be- 
tween £ for the largest spin in our study (S = 5/2) and 
the classical antiferromagnet (S — > oo) is already very 
small. This result supports the hypothesis that the de- 
viations from asymptotic RC behavior evident in Figs.[l] 
and |^ are driven primarily by RC to classical crossover 
effects. 

Another important observable of the scaling theory 
is the Lorentzian amplitude of the spin-spin correla- 
tor, So/£ 2 , where Sq is the correlator magnitude at 
Q = (7r/a,7r/a). This ratio has only a power-law tem- 
perature dependence, and is therefore less sensitive to 
the model parameters than £ or So separately. A de- 
tailed discussion of this quantity is beyond the scope of 
this publication. However, we note some important issues 
which need to be clarified before the scaling behavior of 
So/£ 2 can be fully understood. 

The RC scaling prediction for this quantity is || 
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where Z3 is a universal number and Nq is the T = 
sublattice magnetization defined such that in the classical 
(S — > 00) limit Nq = S, and So is defined such that 
S = S(S + l)/3 for T -> 00. The value of Z 3 can be 
estimated independently using the numerical data for the 
two limiting cases 5=1/2 and S -► 00 

Z3 is estimated by substituting series expansion data 
at the lowest temperature where calculation is possible, 
Tmin, into Eq.(^|). The data so-obtained for S = 1/2 
and S — > 00 differ by a factor of two: Z3 w 3.2 for 
S = 1/2 when estimated at T m i n = 0.35 J, and Z3 « 6.6 
for S — > 00 when estimated at T m in ~ 0.8JS 2 . The 
disagreement of the S — 1/2 and 5 — > 00 results for 
this universal parameter suggests that at least one of the 
models is not in the scaling limit at the respective T m i n 

Here we speculate on two possible causes of this dis- 
crepancy. First, if the classical model is in the classical 
scaling limit and the 5=1/2 model is not in RC limit, 
this discrepancy may be due to the vicinity of the classi- 
cal to RC crossover. In this case, So/£ 2 T 2 for the 5=1/2 
model should increase at temperatures lower than those 
studied numerically, in order to reach its presumed larger 
scaling limit for T — ► 0. This scenario may be consis- 
tent with neutron scattering measurements which give 
S /e ~ const §|H. 

Second, for either of the models the true T — > scaling 
form for large g£ (q is a deviation from Q = (n/a, w/a)), 
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cannot hold simultaneously at the temperatures studied 
for any q inside the Brillouin zone. The reason is that 
the RC scaling form (0), which has to match the hy- 
drodynamic expression (||) for q£ — > 00, is numerically 
found 11 1^] to approach its asymptotic large-q limit 
quite slowly. In the temperature range studied, it still 
has large deviations from this limit for any q inside the 
Brillouin zone (q < a- 1 ), hence Eqs.(@) and (§) are not 
equal and at least one of them must be inconsistent with 
the data. 

However, both Eqs.(Q) and (g) are derived using essen- 
tially the same assumptions that (i) c/£ <C T and that 
(ii) T is much smaller than other energy scales. This 
might indicate the importance of lattice corrections (i.e. 
finite size of the Brillouin zone) to the ratio So /£ 2 in the 
temperature range studied, in which £ is itself consistent 
with our proposed crossover scaling scenario for any 5. 

In summary, we present and analyze high-temperature 
series expansion data for the spin-S square-lattice Heisen- 
berg model for all 5 in the range 5 = 1/2 to 5 = 5/2. 
In agreement with neutron scattering measurements, we 
find that the correlation length deviates from the low- 
temperature RC prediction of Eq.(Q) for 5 > 1/2, and 
that the deviation becomes larger for larger 5. We sug- 
gest that this deviation primarily reflects a crossover to 
classical (as opposed to renormalized classical, or RC) be- 
havior in the temperature range c/a <C T <C p s , or equiv- 
alently JS «C T <C JS 2 , where the correlation length is 
nearly equal to that of the classical 5 — > 00 model. In 
a subsequent publication |]l7f , which is now in prepara- 
tion, we propose to use comparisons between AFMs and 
FMs with the same value of 5 to study the role of quan- 
tum effects as a function of spin; the results obtained in 
]l7t using this method are consistent with the crossover 
scenario proposed here. Finally, we discuss the behav- 
ior of the staggered structure factor Sq, which measures 
the strength of magnetic correlations, and speculate on 
possible causes of its apparent deviation from the T — > 
scaling limit. 
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and the hydrodynamic result, valid for «5«T/c, 
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